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Preface 



The aim of this book is two fold. At the outset the book gives 
most of the available literature about Fuzzy Relational Equations 
(FREs) and its properties for there is no book that solely caters to 
FREs and its applications. Though we have a comprehensive 
bibliography, we do not promise to give all the possible available 
literature about FRE and its applications. We have given only 
those papers which we could access and which interested us 
specially. We have taken those papers which in our opinion could 
be transformed for neutrosophic study. 

The second importance of this book is that for the first time 
we introduce the notion of Neutrosophic Relational Equations 
(NRE) which are analogous structure of FREs. Neutrosophic 
Relational Equations have a role to play for we see that in most of 
the real-world problems, the concept of indeterminacy certainly 
has its say; but the FRE has no power to deal with indeterminacy, 
but this new tool NRE has the capacity to include the notion of 
indeterminacy. So we feel the NREs are better tools than FREs to 
use when the problem under investigation has indeterminates. 
Thus we have defined in this book NREs and just sketched its 
probable applications. 

This book has five chapters. The first chapter is a bulky one 
with 28 sections. These sections deal solely with FREs and their 
properties. By no means do we venture to give any proof for the 
results for this would make our book unwieldy and enormous in 
size. For proofs, one can refer the papers that have been cited in 
the bibliography. 

The second chapter deals with the applications of FRE. This 
has 10 sections: we elaborately give the applications of FRE in 
flow rates in chemical industry problems, preference and 
determination of peak hour in the transportation problems, the 
social problems faced by bonded laborers etc. 

Chapter three for the first time defines several new 
neutrosophic concepts starting from the notion of neutrosophic 
fuzzy set, neutrosophic fuzzy matrix, neutrosophic lattices, 
neutrosophic norms etc. and just indicate some of its important 
analogous properties. This chapter has six sections which are 
solely devoted to the introduction of several neutrosophic 
concepts which are essential for the further study of NRE. 
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Chapter four has eleven sections. This chapter gives all basic 
notions and definitions about the NREs and introduces NREs. 
Section 4.11 is completely devoted to suggest how one can apply 
NREs in the study of real world problems. We suggest many 
problems in chapter five for the reader to solve. 

This is the third book in the Neutrosophics Series. The earlier 
two books are Fuzzy Cognitive Maps and Neutrosophic Cognitive 
Maps ( http://gallup.unm.edu/~smarandache/NCMs.pdf ) and 
Analysis of Social Aspects of Migrant Labourers Living With 
HIV/AIDS Using Fuzzy Theory and Neutrosophic Cognitive 
Maps: With Specific Reference to Rural Tamil Nadu in India 
( http://gallup.unm.edu/~smarandache/NeutrosophvAIDS.pdf ). 

Finally, we thank Meena Kandasamy for the cover design. 
We thank Kama Kandasamy for the layout of the book and for 
drawing all the figures used in this book. She displayed an 
enormous patience that is worthy of praise. We owe deep thanks 
to Dr.K.Kandasamy for his patient proof-reading of the book. 
Without his help this book would not have been possible. 



W.B.VASANTHA KANDASAMY 
FLORENTIN SMARANDACHE 




Chapter One 



FUZZY RELATIONAL EQUATIONS: 
BASIC CONCEPTS AND 
PROPERTIES 



The notion of fuzzy relational equations based upon the max-min 
composition was first investigated by Sanchez [84]. He studied 
conditions and theoretical methods to resolve fuzzy relations on 
fuzzy sets defined as mappings from sets to [0,1]. Some theorems 
for existence and determination of solutions of certain basic fuzzy 
relation equations were given by him. However the solution 
obtained by him is only the greatest element (or the maximum 
solution) derived from the max-min (or min-max) composition of 
fuzzy relations. [84] ’s work has shed some light on this important 
subject. Since then many researchers have been trying to explore 
the problem and develop solution procedures [1, 4, 10-12, 18, 34, 
52,75-80,82,108,111], 

The max-min composition is commonly used when a system 
requires conservative solutions in the sense that the goodness of 
one value cannot compensate the badness of another value [117]. 
In reality there are situations that allow compensatability among 
the values of a solution vector. In such cases the min operator is 
not the best choice for the intersection of fuzzy sets, but max- 
product composition, is preferred since it can yield better or at 
least equivalent result. Before we go into the discussion of these 
Fuzzy Relational Equations (FRE) and its properties it uses and 
applications we just describe them. 

This chapter has 28 sections that deal with the properties of 
FRE, methods of solving FRE using algorithms given by several 
researchers and in some cases methods of neural networks and 
genetic algorithm is used in solving problems. A complete set of 
references is given in the end of the book citing the names of all 
researchers whose research papers have been used. 
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1.1 Binary Fuzzy Relation and their properties 

It is well known fact that binary relations are generalized 
mathematical functions. Contrary to functions from X to Y, binary 
relations R(X, Y) may assign to each element of X two or more 
elements of Y. Some basic operations on functions such as the 
inverse and composition are applicable to binary relations as well. 

Given a fuzzy relation R(X, Y), its domain is a fuzzy set on 
X, dom R, whose membership function is defined by 

dom (R(x)) = max R (x, y) 

yeY 



for each x e X. That is, each element of set X belongs to the 
domain of R to the degree equal to the strength of its strongest 
relation to any member of set Y. The range of R (X, Y) is a fuzzy 
relation on Y, ran R whose membership function is defined by 

ran R(y) = max R (x, y) 

xgX 



for each y e Y. That is, the strength of the strongest relation that 
each element of Y has to an element of X is equal to the degree of 
that elements membership in the range of R. In addition, the 
height of a fuzzy relation R(X,Y) is a number, h(R), defined by 

h (R) = max max (R (x, y). 

ye Y xeX 



That is h(R) is the largest membership grade attained by any pair 
(x, y) in R. 

A convenient representation of binary relation R(X, Y) are 
membership matrices R = [r xy ] where r xy = R(x, y). Another useful 
representation of binary relation is a sagittal diagram. Each of the 
sets X, Y is represented by a set of nodes in the diagram nodes 
corresponding to one set is distinguished from nodes representing 
the other set. 

Elements of X x Y with non-zero membership grades in R(X, 
Y) are represented in the diagram by lines connecting the 
respective nodes. 
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We illustrate the sagittal diagram of a binary fuzzy relation 
R(X, Y) together with the corresponding membership matrix in 
Figure 1.1.1. 




The inverse of a fuzzy relation R(X, Y) denoted by R'*(Y, X) 
is a relation on Y x X defined by R" 1 (y, x) = R (x, y) for all x e X 
and for all yeX. A membership matrix R' 1 = | r 1 yx ] representing 
R" 1 (Y, X) is the transpose of the matrix R for R (X, Y) which 
means that the rows of R" 1 equal the columns of R and the 
columns of R" 1 equal the rows of R. 

Clearly (R 1 )' 1 = R for any binary fuzzy relation. Thus a 
fuzzy binary relation can be represented by the sagittal diagram. 
The corresponding membership matrix: 



yi yi ys y4 ys 



*1 

*2 

*3 

X 4 

*5 

*6 

x 7 



0 

0 

.2 

0 

0 

0 

.2 



.7 .5 0 

.4 0 .1 

0 0 0 
0 .1 1 
0 0 .3 

0 0 .6 
0 .8 0 



0 

0 

0 

0 

.7 

.7 

.5 



R is the membership matrix. 
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Consider now two binary fuzzy relations P(X, Y) and Q(Y, 
Z) with a common set Y. The standard composition of these 
relations, which is denoted by P(X, Y) ° Q(Y, Z), produces a 
binary relation R (X, Z) on X x Z defined by 

R(x,z) = [P ° Q] (x, z) 

max min [P(x, y), Q (y, z)] 

yeY 



for all x e X and all z e Z. This composition, which is based on 
the standard t-norm and t-co-norm is often referred to as the max- 
min composition. It follows directly from the above equation that 

[P (X, Y) ° Q (Y, Z )] 1 = Q " 1 (Z, Y) 0 P 1 (Y, X) 

[P (X, Y) ° Q (Y, Z)] ° R (Z, W) = P (X, Y) 0 [Q (Y, Z) 0 
R (Z, W)]. 

This is the standard (or max-min) composition, which is 
associative, and its inverse is equal to the reverse composition of 
the inverse relations. 

However the standard composition is not commutative 
because Q(Y, Z) ° P(X, Y) is not even well defined when X * Z. 
Even if X = Z and Q (Y, Z) ° P(X, Y) are well defined, we may 
have P(X, Y) ° Q (Y, Z) * Q (Y, Z ) ° P(X, Y). Compositions of 
binary fuzzy relations can be performed conveniently in terms of 
membership matrices of the relations. Let P = [p ik ], Q = [q kj ] and 
R = [ry] be the membership matrices of binary relations such that 
R = P ° Q. We can then write using this matrix notations 

[r,j] = [Pik] ° [qkj] 

where r^ = max min (p ik q kj ). Observe that the same elements of P 

k 

and Q are used in the calculation of R as would be used in the 
regular multiplication of matrices, but the product and sum 
operations are here replaced with max and min operations 
respectively. 

A similar operation on two binary relations, which differs 
from the composition in that it yields triples instead of pairs, is 
known as the relational join. For fuzzy relations P(X, Y) and 
Q(Y, Z), the relational join P * Q, corresponding to the standard 
max-min composition is a ternary relation R(X, Y, Z) defined by 
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R (x, y, z) = [P * Q] (x, y, z) = min [P (x, y), Q (y, z)] for 
each x £ X, y e Y and z £ Z. 

The fact that the relational join produces a ternary relation 
from two binary relations is a major difference from the 
composition, which results in another binary relation. 




Formally [P °Q] (x, z) = max [P * Q] (x, y, z) for each x £ X 
and z £ Z. Now we just see what happens if the binary relation on 
a single set. Binary relation R (X, X) can be expressed by the 
same forms as general binary relations. 

The following properties are to be observed: 

i. Each element of the set X is represented as a single 
node in the diagram. 

ii. Directed connection between nodes indicates pairs 
of elements of X, with the grade of membership in R 
is non-zero. 

iii. Each connection in the diagram is labeled by the 
actual membership grade of the corresponding pair 
in R. 

An example of this diagram for a relation R (X, X) defined on X = 
{xi, x 2 , x 3 , x 4 , x 5 } is shown in Figure 1.1.3. A crisp relation R (X, 
X) is reflexive if and only if (x, x) £ R. for each x £ R, that is if 
every element of X is related to itself, otherwise R(X, X) is called 
irreflexive. If (x, x) g R for every x £ X the relation is called anti 
reflexive. 

A crisp relation R(X, X) is symmetric if and only if for every 
(x, y) £ R, it is also the case that (y, x) £ R where x, y £ X. Thus 
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whenever an element x is related to an element y through a 
symmetric relation, y is also related to x. If this is not the case for 
some x, y then the relation is called asymmetric. If both (x, y) e R 
and (y, x) e R implies x = y then the relation is called anti 
symmetric. If either (x, y) e R or (y, x) e R whenever x^y, then 
the relation is called strictly anti symmetric. 



X l X2 X3 X4 X5 



*1 

*2 

*3 

X 4 

X 5 



.2 

0 

.1 

.2 

.6 



0 .1 0 .7 

.8 .6 0 0 

0 .2 0 0 

0 0 .1 .4 

.1 0 0 .5 



The corresponding sagittal diagram is given in Figure 1.1.3: 




.2 -5 .8 
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Table 



X 


y 


R(x, y) 


Xi 


Xi 


.2 


Xi 


X 3 


.1 


Xi 


X 5 


.7 


X 2 


X 2 


.8 


X 2 


X 3 


.6 


X 3 


Xi 


.1 


X 3 


X 3 


.2 


x 4 


Xl 


.2 


X 4 


X 4 


.1 


X 4 


X 5 


.4 


X 5 


Xi 


.6 


X 5 


X 2 


.1 


X 5 


X 5 


.5 



A crisp relation R (X, Y) is called transitive if and only if (x, z) e 
R, whenever both (x, y) e R and (y, z) e R for at least one y e X. 

In other words the relation of x to y and of y to z imply the 
relation x to z is a transitive relation. A relation that does not 
satisfy this property is called non-transitive. If (x, z) g R 
whenever both (x, y) e R and (y, z) e R, then the relation is 
called anti-transitive. The reflexivity, symmetry and transitivity is 
described by the following Figure 1.1.5: 





Figure: 1 .1 .5 
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A fuzzy relation R (X, X) is reflexive if and only if R(x,x) = 1 
for all x e X, if this is not the case for same x e X, the relation is 
called irreflexive, if it is not satisfied for all x e X, the relation is 
called anti reflexive. A weaker form of reflexivity referred to as 
e - reflexivity denoted by R (x, x) > e where 0 < e < 1. A fuzzy 
relation is symmetric if and only if 

R(x, y) = R (y, x) 

for all x, y e X, if this relation is not true for some x, y e X, the 
relation is called anti symmetric. Further more when R (x, y) > 0 
and R (y, x) > 0 implies x = y for all x, y e X the relation R is 
called anti symmetric. 

A fuzzy relation R (X, X) is transitive if R (x, z) > max min 

y<EY 

[R (x, y), R (y, z)] is satisfied for each pair (x, z) e X 2 . A relation 
failing to satisfy this inequality for some members of X is called 
non-transitive and if 

R (x, z) < max min [R (x, y), R (y, z)] 

yeY 

for all (x,z) e X 2 , then the relation is called anti transitive. 



1.2 Properties of Fuzzy Relations 

In this section we just recollect the properties of fuzzy relations 
like, fuzzy equivalence relation, fuzzy compatibility relations, 
fuzzy ordering relations, fuzzy moiphisms and sup-i-compositions 
of fuzzy relation. For more about these concepts please refer [43]. 

Now we proceed on to define fuzzy equivalence relation. A 
crisp binary relation R(X, X) that is reflexive, symmetric and 
transitive is called an equivalence relation. For each element x in 
X, we can define a crisp set A x , which contains all the elements of 
X that are related to x, by the equivalence relation. 

A x = {y I (x, y) e R (X, X)} 

A x is clearly a subset of X. The element x is itself contained in A x 
due to the reflexivity of R, because R is transitive and symmetric 
each member of A x , is related to all the other members of A x . 
Further no member of A x , is related to any element of X not 
included in A x . This set A x is referred to an as equivalence class 
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of R (X, X) with respect to x. The members of each equivalence 
class can be considered equivalent to each other and only to each 
other under the relation R. The family of all such equivalence 
classes defined by the relation which is usually denoted by X / R, 
forms a partition on X. 

A fuzzy binary relation that is reflexive, symmetric and 
transitive is known as a fuzzy equivalence relation or similarity 
relation. In the rest of this section let us use the latter term. While 
the max-min form of transitivity is assumed, in the following 
discussion on concepts; can be generalized to the alternative 
definition of fuzzy transitivity. 

While an equivalence relation clearly groups elements that 
are equivalent under the relation into disjoint classes, the 
interpretation of a similarity relation can be approached in two 
different ways. First it can be considered to effectively group 
elements into crisp sets whose members are similar to each other 
to some specified degree. Obviously when this degree is equal to 
1, the grouping is an equivalence class. Alternatively however we 
may wish to consider the degree of similarity that the elements of 
X have to some specified element x e X. Thus for each x e X, a 
similarity class can be defined as a fuzzy set in which the 
membership grade of any particular element represents the 
similarity of that element to the element x. If all the elements in 
the class are similar to x to the degree of 1 and similar to all 
elements outside the set to the degree of 0 then the grouping again 
becomes an equivalence class. We know every fuzzy relation R 
can be uniquely represented in terms of its a-cuts by the formula 

R = u a. a R 

ae(0,l] 



It is easily verified that if R is a similarity relation then each a- 
cut, a R is a crisp equivalence relation. Thus we may use any 
similarity relation R and by taking an a - cut a R for any value a e 
(0, 1], create a crisp equivalence relation that represents the 
presence of similarity between the elements to the degree a. Each 
of these equivalence relations form a partition of X. Let it ( a R) 
denote the partition corresponding to the equivalence relation “R. 
Clearly any two elements x and y belong to the same block of this 
partition if and only if R (x, y) > a. Each similarity relation is 
associated with the set ji (R) = {ji ( a R) | a e (0,1]} of partition of 
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X. These partitions are nested in the sense that n ( a R) is a 
refinement of n ( ^R) if and only a > (3. 

The equivalence classes formed by the levels of refinement of 
a similarity relation can be interpreted as grouping elements that 
are similar to each other and only to each other to a degree not 
less than a. 

Just as equivalences classes are defined by an equivalence 
relation, similarity classes are defined by a similarity relation. For 
a given similarity relation R(X, X) the similarity class for each x 
e X is a fuzzy set in which the membership grade of each element 
y e X is simply the strength of that elements relation to x or R(x, 
y). Thus the similarity class for an element x represents the degree 
to which all the other members of X are similar to x. Expect in the 
restricted case of equivalence classes themselves, similarity 
classes are fuzzy and therefore not generally disjoint. 

Similarity relations are conveniently represented by 
membership matrices. Given a similarity relation R, the similarity 
class for each element is defined by the row of the membership 
matrix of R that corresponds to that element. 

Fuzzy equivalence is a cutworthy property of binary relation 
R(X, X) since it is preserved in the classical sense in each a-cut of 
R. This implies that the properties of fuzzy reflexivity, symmetry 
and max-min transitivity are also cutworthy. Binary relations are 
symmetric and transitive but not reflexive are usually referred to 
as quasi equivalence relations. 

The notion of fuzzy equations is associated with the concept 
of compositions of binary relations. The composition of two fuzzy 
binary relations P (X, Y) and Q (Y, Z) can be defined, in general 
in terms of an operation on the membership matrices of P and Q 
that resembles matrix multiplication. This operation involves 
exactly the same combinations of matrix entries as in the regular 
matrix multiplication. However the multiplication and addition 
that are applied to these combinations in the matrix multiplication 
are replaced with other operations, these alternative operations 
represent in each given context the appropriate operations of 
fuzzy set intersections and union respectively. In the max-min 
composition for example, the multiplication and addition are 
replaced with the min and max operations respectively. 

We shall give the notational conventions. Consider three 
fuzzy binary relations P (X, Y), Q (Y, Z) and R (X, Z) which are 
defined on the sets 



X = {X] | i e 1} 
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Y = {yj | j e J} and 
Z = {z k | k e K} 

where we assume that I = N n J = N m and K = N 5 . Let the 
membership matrices of P, Q and R be denoted by P = [pij], Q = 
[q,j], R = M respectively, where p y = P (Xj, yj), q jk = Q (y j5 z k ) r ;j 
= R (x ; , z k ) for all iel (=N n ), jeJ = (N m ) and keK (=N 5 ). This 
clearly implies that all entries in the matrices P, Q, and R are real 
numbers from the unit interval [0, 1], Assume now that the three 
relations constrain each other in such a way that P°Q = R where ° 
denotes max-min composition. This means that max min (p;j qj k ) 

j<^J 

= r ik for all iel and ke" K. That is the matrix equation P° Q = R 
encompasses n x s simultaneous equations of the form 
max min (py, q jk ) = r ik . When two of the components in each of 

j^J 

the equations are given and one is unknown these equations are 
referred to as fuzzy relation equations. 

When matrices P and Q are given the matrix R is to 
determined using P ° Q = R. The problem is trivial. It is solved 
simply by performing the max-min multiplication - like operation 
on P and Q as defined by max min f py, q jk ) = r ik . Clearly the 

jcJ 

solution in this case exists and is unique. The problem becomes 
far from trivial when one of the two matrices on the left hand side 
of P ° Q = R is unknown. In this case the solution is guaranteed 
neither to exist nor to be unique. 

Since R in P ° Q = R is obtained by composing P and Q it is 
suggestive to view the problem of determining P (or alternatively 
Q ) from R to Q (or alternatively R and P) as a decomposition of 
R with respect to Q (or alternatively with respect to P). Since 
many problems in various contexts can be formulated as problems 
of decomposition, the utility of any method for solving P°Q = R 
is quite high. The use of fuzzy relation equations in some 
applications is illustrated. Assume that we have a method for 
solving P°Q = R only for the first decomposition problem (given 
Q and R). 

Then we can directly utilize this method for solving the 
second decomposition problem as well. We simply write P ° Q = 
R in the form Q' 1 o P" 1 = R" 1 employing transposed matrices. We 
can solve Q" 1 o P" 1 = R" 1 for Q" 1 by our method and then obtain the 
solution of P ° Q = R by (Q 1 )" 1 = Q. 
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We study the problem of partitioning the equations P ° Q = R. 
We assume that a specific pair of matrices R and Q in the 
equations P ° Q = R is given. Let each particular matrix P that 
satisfies P ° Q = R is called its solution and let S (Q, R) = (P | P ° 
Q = R} denote the set of all solutions (the solution set). 

It is easy to see this problem can be partitioned, without loss 
of generality into a set of simpler problems expressed by the 
matrix equations p; o Q = q for all iel where 
p i = [Pij I j e J] and 
q = [r ik | k e K], 

Indeed each of the equation in max min (pijqjk) = q k contains 

j<^J 

unknown p ;j identified only by one particular value of the index i, 
that is, the unknown p,j distinguished by different values of i do 
not appear together in any of the individual equations. Observe 
that pi, Q, and q in p ; ° Q = q represent respectively, a fuzzy set on 
Y, a fuzzy relation on Y x Z and a fuzzy set on Z. Let S; (Q, q) = 
[Pi I Pi o Q = q] denote, for each iel, the solution set of one of the 
simpler problem expressed by pi ° Q = q. 

Thus the matrices P in S (Q, R) = [P | P ° Q = R ] can be 
viewed as one column matrix 



Pi 



P = 



P 2 



P 



n 



where pi e S; (Q, q) for all i e I = (=N n ). It follows immediately 
from max min (pij q jk ) = r ik . That if max q jk < r ik for some i e I 

j*J i* J 

and some k e K, then no values p;j e [0, 1] exists (j e J) that 
satisfy P ° Q = R, therefore no matrix P exists that satisfies the 
matrix equation. This proposition can be stated more concisely as 
follows if 



max 

i^J 



Vjk 






for some k e K then S (Q, R) = 4>. This proposition allows us in 
certain cases to determine quickly that P°Q = R has no solutions 
its negation however is only a necessary not sufficient condition 
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for the existence of a solution of P ° Q = R that is for S (Q, R) * 4». 
Since P°Q = R can be partitioned without loss of generality into 
a set of equations of the form p; ° Q = q we need only methods for 
solving equations of the later form in order to arrive at a solution. 
We may therefore restrict our further discussion of matrix 
equations of the form P ° Q = R to matrix equation of the simpler 
form P ° Q = r, where p = [pj | j e J], Q = [q,k | j e J, k 6 K] and 
r = (r k | k e K], 

We just recall the solution method as discussed by [43]. For 
the sake of consistency with our previous discussion, let us again 
assume that p, Q and r represent respectively a fuzzy set on Y, a 
fuzzy relation on Y x Z and a fuzzy set on Z. Moreover let J = N m 
and K = N s and let S (Q, r) = {p | p ° Q = r} denote the solution set 
of 

P°Q = r. 

In order to describe a method of solving p ° Q = r we need to 
introduce some additional concepts and convenient notation. First 
let p denote the set of all possible vectors. 

P = {Pj I j e J} 

such that pj e [0, 1] for all j € J, and let a partial ordering on 
p be defined as follows for any pair p 1 , p 2 e p p 1 < p 2 if and 
only if p 2 < pj for all j eJ. Given an arbitrary pair p 1 , p 2 e p 

such that p 1 < p 2 let [p 1 ,p 2 ]={pe p | p 1 < p < p 2 }. For any pair 
p 1 , p 2 e p ({p 1 , p 2 } < } is a lattice. 

Now we recall some of the properties of the solution set S (Q, 
r). Employing the partial ordering on p, let an element p of S (Q, 
r) be called a maximal solution of p ° Q = r if for all p e S (Q, r), 
p > p implies p ■*= p if for all p e S (Q, r) p < p then that is the 
maximum solution. Similar discussion can be made on the 
minimal solution of p 0 Q = r. The minimal solution is unique if 
p > p (i.e. p is unique). 

It is well known when ever the solution set S (Q, r) is not 
empty it always contains a unique maximum solution p and it 

may contain several minimal solution. Let S (Q, r) denote the set 
of all minimal solutions. It is known that the solution set S (Q, r) 
is fully characterized by the maximum and minimal solution in 
the sense that it consists exactly of the maximum solution p all 
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the minimal solutions and all elements of p that are between p 
and the numeral solution. 

Thus S (Q, r) — u [p, p\ where the union is taken for all 

p 

p gS (Q, r). When S (Q, r) * (|>, the maximum solution. 

P ~ t P j I j e J] of p ° Q = r is determined as follows: 

P j = min ct (q ik , r k ) where ct (q jk , r k ) = 

J ksK 

when p determined in this way does not satisfy p ° Q = r then 
S(Q, r) = <j). That is the existence of the maximum solution p as 
determined by p , = min ct (q; k , r k ) is a necessary and sufficient 

1 keK 

condition for S (Q, r) * <b. Once p is determined by p , = min ct 

1 ksK 

(q; k , r k ), we must check to see if it satisfies the given matrix 
equations p ° Q = r. If it does not then the equation has no solution 
(S (Q, r) = 4>), otherwise p in the maximum solution of the 

equation and we next determine the set S (Q, r) of its minimal 
solutions. 



r k if Pjk > r k 
1 otherwise 



1.3 Fuzzy compatibility relations and composition of fuzzy relations 

In this section we recall the definition of fuzzy compatibility 
relations, fuzzy ordering relations, fuzzy morphisms, and sup and 
inf compositions of fuzzy relations. 

DEFINITION [43]: A binary relation R(X, X) that is reflexive and 
symmetric is usually called a compatibility relation or tolerance 
relation. When R(X, X) is a reflexive and symmetric fuzzy relation 
it is sometimes called proximity relation. 

An important concept associated with compatibility relations is 
compatibility classes. Given a crisp compatibility relation R(X, 
X), a compatibility class is a subset A of X such that (x, y) e R 
for all x, y e A. A maximal compatibility class or maximal 
compatible is a compatibility class that is not properly contained 
with in any other compatibility class. The family consisting of all 
the maximal compatibles induced by R on X is called a complete 
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cover of X with respect to R. When R is a fuzzy compatibility 
relation, compatibility classes are defined in terms of a specified 
membership degree a. An a-compatibility class is a subset A of X 
such that R (x, y) > a for all x, y e A. Maximal a-compatibles 
and complete a-cover are obvious generalizations of the 
corresponding concepts for crisp compatibility relations. 

Compatibility relations are often conveniently viewed as 
reflexive undirected graphs contrary to fuzzy cognitive maps that 
are directed graphs. In this context, reflexivity implies that each 
node of the graph has a loop connecting the node to itself the 
loops are usually omitted from the visual representations of the 
graph although they are assumed to be present. Connections 
between nodes as defined by the relation, are not directed, since 
the property of symmetry guarantees that all existing connections 
appear in both directions. Each connection is labeled with the 
value corresponding to the membership grade R (x, y) = R (y,x). 

We illustrate this by the following example. 

Example 1.3.1: Consider a fuzzy relation R(X, X) defined on X 
= {X|, x 2 ,..., x 8 } by the following membership matrix: 



X i Xy Xj X4 X$ Xfr Xy Xg 



X 1 
*2 
*3 
X 4 

*5 

*6 

*7 

*8 



1 

.3 

0 

0 

.4 

0 

0 

.6 



.3 

1 

.5 

.3 

0 

0 

0 

0 



0 0 
.5 .3 
1 0 
0 1 
0 .2 
.7 0 
.6 .7 
.8 .5 



.4 

0 

0 

.2 

1 

0 

0 

0 



0 0 .6 

0 0 0 

.7 .6 .8 

0 .7 .5 

0 0 0 

1 .2 0 

.2 1 .8 

0 .8 1 



Since the matrix is symmetric and all entries on the main diagonal 
are equal to 1, the relation represented is reflexive, and symmetric 
therefore it is a compatibility relation. The graph of the relation is 
shown by the following figure 1.3.1, its complete a-covers for a > 
0 and a e {0, .3, .1, .4, .6, .5, .2, .8, .7, 1} is depicted. Figure 
1.3.1 is the graph of compatibility relation given in example 1.3.1 
while similarity and compatibility relations are characterized by 
symmetry, ordering relations require asymmetry (or anti- 
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symmetry) and transitivity. There are several types of ordering 
relations. 




A crisp binary relation R(X, X) that is reflexive, anti 
symmetric and transitive is called a partial ordering. The common 
symbol < is suggestive of the properties of this class of relations. 
Thus x < y denotes (x, y) e R and signifies that x precedes y. The 
inverse partial ordering R' 1 (X, X) is suggested by the symbol >. 

If y > x including that (y, x) e R' 1 then we say that y 
succeeds x. When x < y; x is also referred to as a predecessor of y 
while y is called a successor of x. When x < y and there is no z 
such that x < y and z < y, x is called an immediate predecessor of 
y and y is called an immediate successor of x. If we need to 
distinguish several partial orderings, such as P, Q and R we use 

p Q R 

the symbol < , < and < respectively. 

Observe that a partial ordering ‘<’ on X does not guarantee 
that all pairs of elements x, y in X are comparable in the sense that 
either x < y or y < x. Thus, for some x, y e X it is possible that x 
is neither a predecessor nor a successor of y. Such pairs are called 
non comparable with respect to <. 
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The following are definitions of some fundamental concepts 
associated with partial orderings: 

1. Ifx e X and x < y for every y e X then x is called the 
first member (or minimum) of X with respect to the 
relation denoted by <. 

2. Ifx s X and y < x for every y e X, then x is called the 
last member (or maximum) of X with respect to the 
partial ordering relation. 

3. If x e X and y < x implies x = y, then x is called a 
minimal member of X with respect to the relation. 

4. If x £ X and x < y implies x = y, then x is called a 
maximal member of X with respect to the relation [43]. 
Using these concepts every partial ordering satisfies the 
following properties: 

i. There exist at most one first member and at 
most one last member. 

ii. There may exist several maximal members and 
several minimal members. 

iii. If a first member exists then only one minimal 
member exists and it is identical with the first 
member. 

iv. If a last member exists, then only one maximal 
member exists and it is identical with the last 
member. 

v. The first and last members of a partial ordering 
relation correspond to the last and first members 
of the inverse partial ordering, respectively. 

Let X again be a set on which a partial ordering is defined and let 
A be a subset ofX(Ac X). Ifx £ X and x < y for every y e A, 
then x is called a lower bound of A on X with respect to the 
partial ordering. Ifx £ X and y < x for every y e A, then x is 
called an upper bound of A on X with respect to the relation. If a 
particular lower bound succeeds, every other lower bound of A, 
then it is called the greatest lower bound or infimum, of A. If a 
particular upper bound proceeds every other upper bound of A 
then it is called the least upper bound or superemum of A. 
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A partial ordering on a set X that contains a greatest lower 
bound and a least upper bound for every subset of two elements of 
X is called a lattice. 

A partial ordering < on X is said to be connected if and only 
if for all x, y e X (x * y) implies either x < y or y < x. When a 
partial ordering is connected all pairs of elements of X are 
comparable by the ordering such an ordering is usually called a 
linear ordering, some alternative names used in the literature are 
total ordering, simple ordering and complete ordering. 

Partial ordering can be represented by diagrams and this sort 
of diagrams are called Hasse’s diagrams. 

A fuzzy binary relation R on a set X is a fuzzy partial 
ordering if and only if it is reflexive anti-symmetric and transitive 
under some form of fuzzy transitivity. Any fuzzy partial ordering 
based on max-min transitivity can be resolved into a series of 
crisp partial orderings in the same way in which this is done for 
similarity relations, that is by taking a series of a-cuts that 
produce increasing levels of refinement. When a fuzzy partial 
ordering is defined on a set X, two fuzzy sets are associated with 
each element x in X. The first is called the dominating class of x. 

It is denoted by 

R> [ x ] and is defined by R >[ x] (y) = R (x, y) 

where y e X. In other words the dominating class of x contains 
the members of X to the degree to which they dominate x. The 
second fuzzy set of concern is the class dominated by x, which is 
denoted by 

R< [ x] and defined by R< w (y) = R (y, x) 
where yeX. The class dominated by x contains the elements of 
X to the degree to which they are dominated by x. An element x e 
X is undominated if and only if R (x, y) = 0 for all y e X and x * 
y, an element x is undominating if and only if R (y,x) = 0 for all y 
e X and y * x. For a crisp subset A of a set X on which a fuzzy 
partial ordering R is defined, the fuzzy upper bound for A is the 
frizzy set denoted by U (R, A) and defined by 

U (R, A) = f) %v] 

xeA 

where n denotes an appropriate fuzzy intersection. This definition 
reduces to that of the conventional upper bound when the partial 
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ordering is crisp. If a least upper bound of the set A exists it is the 
unique element x in U (R, A) such that 

U (R, A) (x) > 0 and 

R (x, y) > 0 for all elements y in the support of U (R, A). Several 
other concepts of crisp orderings easily generalize to the fuzzy 
case. A fuzzy preordering is a fuzzy relation that is reflexive and 
transitive. Unlike a partial ordering, the preordering is not 
necessarily anti-symmetric. 

A fuzzy weak ordering R is an ordering satisfying all the 
properties of a fuzzy linear ordering except anti-symmetry. 
Alternatively it can be thought of as a fuzzy preordering in which 
either R (x, y) > 0 or R (y, x) > 0 for all x * y. A fuzzy strict 
ordering is anti-reflexive anti-symmetric and transitive, it can 
clearly be derived from any partial ordering R by replacing the 
values R(x, x) = 1 with zeros for all x. 

Now we proceed on to recall the definition of fuzzy 
morphisms. 

If two crisp binary relations R (X, X) and Q (Y, Y) are 
defined on sets X and Y, respectively then a function h : X — > Y is 
said to be a homomorphism from (X, R) to (Y, Q) if (xj, x 2 ) e R 
implies (h(xj), h(x 2 )) e Q for all x b x 2 e X. In other words a 
homomorphism implies that for every two elements of the set X 
which are related under the relation R, their homomorphic images 
h (xj), h(x 2 ) in set Y are related under the relation Q. 

When R (X, X) and Q (Y, Y) are fuzzy binary relations this 
implication can be generalized to R (xi, x 2 ) < Q (h(xi), h(x 2 )), for 
all x b x 2 e X and their images h (xj), h (x 2 ) e Y. Thus, the 
strength of relation between two elements under R is equated or 
exceeded by the strength of relation between their homomoiphic 
images under Q. 

Note that it is possible for a relation to exist under Q between 
the homomorphic images of two elements that are themselves 
unrelated under R. When this is never the case under a 
homomorphic function h, the function is called a strong 
homomorphism. It satisfies the two implications 

(x b x 2 ) e R implies (h(xi), h(x 2 )> e Q 

for all xi, x 2 e X and (yi, y 2 ) e Q implies (x 1? x 2 ) e R, for all yi, 
y 2 e Y where Xj e h" 1 (yi) and x 2 eh' 1 (y 2 ). Observe that when h 
is many to -one the inverse of h for each element of Y is a set of 
elements from X instead of a single element of X. If relations R 
(X, X) and Q (Y, Y) are fuzzy, then the criteria that a many-to- 
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one function h must satisfy in order to be a strong homomorphism 
are somewhat modified. The function h imposes a partition. 7i h on 
the set X such that any two elements x b x 2 e X belong to the 
same block of the partition if and only if h maps them to the same 
element of Y. Let A = {ai, a 2 , ..., a n } and B = {bi, b 2 ,..., b n } be 
two blocks of this partition 7Xh and let all elements of A be mapped 
to some element yi e Y and all elements of B to some element y 2 
e Y. Then the function h is said to be a strong homomorphism 
from (X, R) to (Y, Q) if and only if the degree of the strongest 
relation between any element of A and any element of B in the 
fuzzy relation R equals the strength of the relation between yi and 
y 2 in the fuzzy relation Q. Formally 

max R(a i b j ) = Q(y,,y 2 ). 

y 

This equality must be satisfied for each pair of blocks of the 
partition ji h . If a homomorphism exists from (X, R) to (Y, Q) 
then the relation Q (X, Y) preserves some of the properties of the 
relation R (X, X) - namely, that all the pairs (x b x 2 ) e X x X 
which are members of R have corresponding homomorphic 
images (h^), h(x 2 )) e Y x Y which are members of Q. Other 
members of Q may exist, however that are not the counterparts of 
any number of R. This is not the case when the homomorphism is 
strong. Here more properties of the relation R are preserved in 
relation Q. In fact Q represents a simplification of R in which 
elements belong to the same block of the partition TCh created by 
the function h on the set X are no longer distinguished. These 
functions are useful for performing various kinds of 
simplifications of systems that preserve desirable properties in 
sets such as ordering or similarity. 

If h : X — > Y is a homomorphism from (X, R) to (Y, Q) and if 
h is completely specified, one to one and on to then it is called as 
isomorphism. This is effectively a translation or direct relabeling 
of elements of the set X into elements of the set Y that preserves 
all the properties of R in Q. If Y c X then h is called an 
endomorphism. A function that is both an isomorphism and an 
endomorphism is called an automorphism. In this case the 
function maps the set X to itself and the relation R and Q are 
equal. Each of these terms applies without modification to fuzzy 
homomorphisms. 

Now we just recall the notions of Sup-i-compositions of 
fuzzy relations. Sup-i compositions of binary fuzzy relations 
where i refers to a t-norm generalizes the standard max-min 
composition. The need to study these generalizations concepts 
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from some applications such as approximate reasoning and fuzzy 
control. Given a particular t-norm i and two fuzzy relations P (X, 
Y) and Q (Y, Z), the Sup-i-composition of P and Q is a fuzzy 

relation P 1 ° Q on X, Y, Z defined by 

(P 1 ° Q} (x, z) = Sup i [P(x, y), Q(y,z)] 

y<=Y 

for all x e X, z e Z. When i is chosen to be the min operator P l a Q 

becomes the standard composition P°Q. 

Given fuzzy relations P(X, Y), P,(X, Y), Q(Y, Z), Qj(Y, Z) 
and R(Z, V) where j takes values in an index set J, the following 
are basic properties of the Sup-i composition under the standard 
fuzzy union and intersections. 

(pi q) i r=p i (q i r) 

f \ 

?! U Q j (p iQ ) 

f \ 

p ‘ ne, e n< p ‘2j) 

\JPjiQ = {J(Pji o Q ) 

j^J i^J 

f)pji o Q s ^iQ ) 

j^j j&j 

(P i Q) _1 = Q 1 i P' 1 

o o 

These properties follow directly from the corresponding 
properties of t-norms and their verification is left for the reader. 

Sup i composition is also monotonic increasing that is for any 
fuzzy relation P(X, Y), Q, (Y, Z), Q 2 (Y, Z) and R(Z, V) if Qi c 
Q then 

P i Qi cP i Q 2 

O O 

Ch i RcQ, i R. 

O O 

The set of all binary fuzzy relations on X 2 forms a complete 
lattice ordered monoid (3(X 2 ), n, u i ) where n and u represent 
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the meet and join of the lattice respectively and i represents, the 

O 

semi group operation of the monoid. The identity of i is defined 

O 

by the relations 

1 1 when x = y 
E (x, y) = \ 

[0 when x ^ y 

The concept of transitivity of a fuzzy relation, which is 
introduced in terms of the max-min composition can be 
generalized in terms of the Sup-i-compositions for the various t- 
norms i. We say the relation R on X 2 is i-transitive if and only if 

R(x, z) > i[R(x, y), R(y, z)J 

for all x, y, z e X. It is easy to show that a fuzzy relation R on X 2 
is i-transitive if and only if R i'RcR which may be used as an 

O 

alternative definition of i-transitivity. 

When a relation R is not i-transitive, we define its i-transitive 
closure as a relation R x <i) that is the smallest i-transitive relation, 
containing R. To investigate properties of i-transitive closure let 

R (n) — R i R (n ~ n 

O 

n = 2, 3,... where R is a fuzzy relation on X 2 and R (l1 = R. Using 
this notation the reader is expected to prove the following 
theorem: 

THEOREM [43]: For any fuzzy relation R on X 2 , the fuzzy relation 

00 

R x (i> = [J R l ' n is the i-transitive closure of R. 

n = 1 

Prove if R be a reflexive fuzzy relation on X 2 , where | X | = n > 2 
then R x (i) = R <n_1) . 

Now we proceed on to recall the notion of inf-Wi 
compositions of fuzzy relations. Give a continuous t-norm i, let 
w ; (a, b) = sup {x e [0, 1] / i(a, x) < b} for every a, b £ [ 0, 1]. 
This operation referred to as operation w, plays an important role 
in fuzzy relation equations. While the t-norm i may be interpreted 
as logical conjunction, the corresponding operation W; may be 
interpreted as logical implication. The following basic properties 
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of W; are left as an exercise for the reader to prove. For any a, aj, 
b, d e [0, 1] where j takes values from an index set J prove W; has 
the following properties: 



1. i (a, b)<d iff Wj(a, d)>b 

2. Wj (w; (a, b) > a 

3 . Wj (i (a, b), d) = w ; (a, w ; (b, d)) 

4. a < d implies 

W; (a, d) > W; (b, d) and 
w; (d, a) < w; (d, b) 

5. i [Wj (a, b), w; (b, d)] < Wj (a, d) 

6. W; [inf aj, b] > sup Wj (aj, b) 

i^J 

7. W; (sup aj, b) = inf W; (a,, b) 

j^J 

8. Wj [b, Sup aj] > sup w; (b, aj) 

j^J j<Ej 

9. Wj (b, inf a f ) = inf w; (b, aj) 

jcJ j^J 

10. i [a, W; (a, b)] < b. 

Prove for the fuzzy relations P(X, Y), Q (Y, Z), R(X, Z) and S(Z, 
V) the following are equivalent 

Pi QcR 

O 

QcP' 1 w, R 

O 

Pc(Q w ; R- 1 ) 1 

O 

Prove P w, (Q w, S) = (P w i Q) w ; S. 

o o o o 

Let P (X, Y), Pj (X, Y), Q(Y, Z) and Qj (Y, Z) be fuzzy 
relations where j takes values in an index set J then prove. 



f ) ( 

U /J , 

\jeJ J o 

f \ 

f] p j 

j o je J o 

f \ 

p w i n Qj =r\ pw ‘ Q j 

° \ jeJ 7 j^J ° 
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p*i\jQj^U Pw iQj- 

° i^J jcJ ° 



Let P (X, Y), Qi(Y, Z), Q 2 (Y, Z) and R (Z, V) be fuzzy relations. 
If Qi £ Q 2 then prove. 

PWjQicP Wj Q 2 and 

O O 

Qi Wj RdQ 2 Wj R. 

o o 

Now if P (X, Y), Q (Y, Z) and R (X, Z) be fuzzy relations prove 
P' 1 i (P Wj Q)cQ 

0 o 

RcPw,. (P _1 i R) 

o 0 

Pc(P Wj Q) vv, Q 1 

O O 

Rc(R Wj Q" 1 ) Wj Q. 



1.4 Optimization Of FRE with Max-Product Composition 

Jiranut Loelamonphing and Shu Cheng Fang [58] in their paper 
“optimization of fuzzy relation equation with max-product 
composition” (2001) has studied the solution set of fuzzy relation 
equations with max product composition and an optimization 
problem with a linear objective function subject to such FRE. By 
identifying the special properties of the feasible domain they 
determine an optimal solution without explicitly generating the 
whole set of minimal solutions. 

The notion of FRE based upon the max-min composition was 
first investigated by [84]. Fie studied conditions and theoretical 
methods to resolve fuzzy relations on fuzzy sets defined as 
mappings from sets into complete Brouwerian lattices. Some 
theorems for existence and determination of solutions of certain 
basic fuzzy relation equations were presented in his work. 
Flowever, the solution obtained in that work is only the greatest 
element (or the maximum solution) derived from the max-min (or 
min-max) composition of fuzzy relations. 

Sanchez’s work [84] has shed some light on this important 
subject. Since then, researchers have been trying to explore the 
problem and develop solution procedures [1, 11, 16, 24, 30, 82]. 
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The “max-min” composition [117] is commonly used when a 
system requires conservative solutions in the sense that the 
goodness of one value cannot compensate the badness of another 
value. In reality, there are situations that allow compensatability 
among the values of a solution vector. In this case, the min 
operator is not the best choice for the intersection of fuzzy sets. 
Instead, the “max-product” composition is preferred since it can 
yield better, or at least equivalent, results. Note that when the 
intersection connector acts non-interactively, it can be uniquely 
defined by the min connector, but when the connector is 
interactive, it is application dependent and cannot be defined 
universally. Some outlines for selecting an appropriate connector 
has been provided by [112, 113]. 

Recently, researchers extended the study of an inverse 
solution of a system of FRE with max-product composition. They 
provided theoretical results for determining the complete solution 
sets as well as the conditions for the existence of resolutions. 
Their results showed that such complete solution sets can be 
characterized by one maximum solution and a number of minimal 
solutions. Since the total number of minimal solutions has a 
combinatorial nature in terms of the problem size, an efficient 
solution procedure is always in demand. 

Motivated by the work of [24], we are interested in studying 
the optimization problem with a linear objective function subject 
to a system of FRE with the “max-product” composition. 

Let A = [ay], 0 < ay < 1, be an (m x n) - dimensional fuzzy 
matrix, b = (bj, ..., b n ), 0 < bj < 1, be an n-dimensional vector, I = 
{1,2, ..., m} and J = {1,2, ..., n}. A system of FRE defined by A 
and b is denoted by 

X o A = b, (1) 

where “o” represents the max-product composition. The 
resolution of (1) is a set of solution vectors x = (x b ..., x m ), 0 < Xj 
< 1, such that 

max {x ; . ay} = bj for j e J (2) 



Let c = (ci,..., c m ) e R m be an m-dimensional vector where C; 
represents the weight (or cost) associated with variable x b for 
i £ I. The optimization problem we are interested in has the 
following form: 
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Minimize 



Z= Yj C ‘ X i (3) 

1=1 



Subject to x o A = b, 

0 < Xj < 1 . 

Note that the characteristics of the solution sets obtained by 
using the max-min operator and the max-product operator are 
similar, i.e., when the solution set is not empty, it can be 
completely determined by a unique maximum solution and a finite 
number of minimal solutions [11, 34]. Since the solution set can 
be non-convex, traditional linear programming methods, such as 
the simplex and interior-point algorithms, cannot be applied to 
this problem. 

[58] denote the solution set of problem (1) by X (A, b) = 
{(x l5 ..., x m ) | x ; e [0, 1], i e I, and x o A = b}. Define X = [0, l] m . 
For x 1 , x 2 e X, we say x 1 < x 2 if and only if x[ < x ■ , Vi e I . In 

this way, the operator “<” forms a partial order relation on X and 
(X, <) becomes a lattice. 

x e X (A, b) is called a maximum solution if x < X for 
all x £ X (A, b). Also x e X (A, b) is called minimal solution if 
x < Jc , for any x e X (A, b) implies x = x . When X (A, b) is not 
empty, it can be completely determined by a unique maximum 
solution and a finite number of minimal solutions [11, 34]. 

The maximum solution can be obtained by applying the 
following operation [58]: 



where 



= A®b = 



A (aii 



0bj) 



iel 



fly 0 bj = 



{ b J /a ij 



'f “a < h i 
V a u >b i 



a a b = min (a, b). 



(4) 



(5) 



Denote the set of all minimal solutions by X (A, b), the complete 
set of solution, X (A, b), is obtained by 
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( 6 ) 



X (A, b) = {x e X | x < x < a} 

x<eX(A, b) 

DEFINITION 1.4.1: For a solution x eX (A, b), we call a 
binding variable if X jo . Cl i ■ = bjfor i 0 £ I and X, a,j < bj, for all i 
£ I. 

When the solution set of (1) is not empty, i.e., X (A, b) &</), we 
define 

I, = { i £l . I a (j = bf, Vj £j, (7) 

A = f xl 2 x...xl n . (8) 

Here /, corresponds to a set of Xj’s that can satisfy constraint j of 
the fuzzy relation equations. And, the set A represents all 
combinations of the binding variables such that every 
combination can satisfy every fuzzy relation constraint. Let each 
combination be represented by f = (f h f 2 , ..., f„) £ A, with f £ f, 
Vj£J. 

The optimization problem (3) can be decomposed into two 
problems, namely 

m 

Minimize Z^^c'jX; (9) 

i=i 

Subject to x o A = b, 

0 < X; < 1 

and 

m 

Minimize Z"=^c" i x i (10) 

i=i 

Subject to x o A = b, 

0 < x; < 1 

where c’ = (c’i, c’ 2 ,..., c’ m ) and c” = (c”i, c” 2 ,..., c” m ) are defined 
such that , Vi el. 
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C, if Ci > 0, 

0 if Ci < 0, 



( 11 ) 



0 if Ci > 0, 

Ci if Ci < 0. 



Apparently, the cost vector c = c’ + c” and the objective value Z’ 
— Z’ + Z”. Intuitively, when all the costs are non-positive, since 
Xj’s are non-negative and the problem is to minimize the objective 
value, we should make x, as large as possible. 

Taking advantage of the special structure studied in the 
previous section, we now introduce some procedures to reduce the 
size of the original problem so that the effort to solve the problem 
is minimized. The key idea behind these reduction procedures is 
that some of the xfs can be determined immediately without 
solving the problem but just by identifying the special 
characteristic of the problem. Special cases which we can 
eliminate from consideration are as follows: 



Case I: q < 0. 

We know that x*j = x t , if C; < 0. Hence, we can take any part that 
are associated with these x t ’s out of consideration. 

Here we define: 



II 

m 

_p 

IA 

© 


(12) 


J = {j e J I Xj, ay =bj, V i e / }. 


(13) 



In other words, / is a set of indices of constraints which can 
be satisfied by a set of x, ’s for i el. Having defined J and / , 
we now eliminate row i, i el and j, j e / from matrix A as well 

as the j th element j £ J , from vector b. 

Let A’ and b’ be the updated fuzzy matrix and fuzzy vector, 
respectively. 

Define J’ = J \ ./ , T represents a reduced set of constraints. 
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Case II : Ij has only one element. 

Consider constraint j e J’. If Ij contains only one element, it 
means that only one x B i e Ij, can satisfy the jth constraint. We 
have x; = fy/ay. 

Define 



1 = {islj | |/j =1; 7'e/'} 


(14) 


J = {jeJ'\x i A i j =b ] ; i G 7} . 


(15) 



Again, we can eliminate row i, i e I , and column j, j e J , 
from the updated fuzzy matrix A’ as well as the j th element j e J , 
from the updated vector b’. Let A” and b” be the reduced fuzzy 
matrix and fizzy vector corresponding to A’ and b’ respectively. 
We also need to update A. Define J”= J’\ J , J” is an index set of 
constraints which need to be solved later by the branch-and-bound 
(B&B) method. The updated A” = IIj e j» Ij. 

The branch-and-bounded algorithm will be performed on 
these A” and b”. If b” is empty, then all constraints have been 
taken care of. Therefore, in order to minimize the objective value, 
since we are now left with positive Cj’s, we can assign the 
minimum value, i.e. zero, to all X;’s whose values have not been 
assigned yet. When b” is not empty, we need to proceed further. 
Details will be discussed in the following: 

In order to identity whether the problem is decomposable, 
consider a set of constraints, say B, which can be satisfied by a 
certain set of variables, say X B . If the decision to choose which 
variable in the set X B to satisfy a constraint in B does not impact 
the decision on the rest of the problem, then we can extract this 
part from the whole problem. 

Let k be the number of sub-problems, 1 < k < ||/"|| . 

Define 



Q 


= {Ijlj 

r 


£ J”}, 

j 


(16) 




4 


rv, 


(17) 




i 


j^j" j 




Q, 


n Q r 




(18) 


Q 


= 0]U 


Q 2 u...u Q k , 


(19) 
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( 20 ) 



a >= n^. 

I (1) = { i | ie Ij, Ij eQ[}, (21) 

J (1) = {j | Ij eOi}. (22) 

In this way, Qi contains sets of Ij’s which have some 
element(s) in common and we can decompose the original 
problem into k sub-problems. 1 (1) and J 1 1 1 correspond to sets of 
indices of variables and constraints, respectively, on which the 
B&B method is performed for sub-problem 1 [58]. 

Problem (9) can be transformed into the following 0-1 integer- 
programming problem 



minimize 



subject to 



z'=Z 

i=l 
m 

Z x ij = 1 



/ 1 


[ bl 1 


\ 


c'; max<^ 


’ X ij 

i a u J 




jeJ 




v 1 


/ 


= 1 


VjeJ 





Xij = 0 or 1 V i e I, j e J. 
Xy = 0 V i, j with i e j. 



(23) 



Algorithm 1 (Algorithm for finding an optimal 
solution) 

Step 1 : Find the maximum solution of ( 1). 

Compute x = A® b = [a" = 1 (a,y 0 bj )] (e/ J according to (4). 



Step 2: Check feasibility. 

If x o A = b, continue. Otherwise, stop! X (A, b) = <)> and problem 
(3) has no feasible solution. 

Step 3: Compute index sets. 

Compute Ij = | i e 1 1 x i ,< 2 y - = bj j, V j e J, which represents a set 
of Xj’s that can satisfy constraint j of the FRE. 
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Step 4: Arrange cost vector. 

Define c’ and c’ according to (11). 



Step 5: Perform problem reduction. 

Compute 7 = {i g I \ c t < 0} and J = {j g J | Xj.a^ = by, i g 7} . 

Eliminate row i e 7 , and column j, j e J , from matrix A to 
obtain A’. Also eliminate the jth element, j e J , from vector b to 
obtain b’. Assign an optimal value x* ; = x, , for i e 7 . If b’ is 
empty, assign zero to unassigned x ; = x t , for i e 7 . If b’ is empty 
assign zero to unassigned x* ; and go to Step 11. Otherwise, 
compute J’ = J\ J and proceed to the next step. 

Step 6: Find singleton Ij. 

Compute 7 = {i e / ; - 1 1|7 ; || = 1; j g /'} and 7 = (j £ J’ | x ; . ay = 

bj; i £ 7 }. Eliminate row i, i e 7 and column j, j e / from 
matrix A’ to obtain A”. Also, eliminate the jth element, j e J , 
from vector b’ to obtain b”. Assign x*j — 'bj / ay, for i e 7 and i e 
Ij. If b” is empty assign zero to unassigned x : and go to Step 11. 
Otherwise, compute J” = J’\ 7 and proceed to the next step. 

Step 7: Decompose the problem. 

Decompose the problem by computing equations (16-22). 

Step 8: Define sub-problems. 

For each sub-problem 1, define problem (11) and its 
corresponding 0-1 inter program using (23). 

Step 9: Solve the integer program(s). 

Solve each integer program by using the branch-and-bound 
method. 

Step 10: Generate an optimal solution of sub-problem. 

For each sub-problem 1, define f 1 = (fj), j £ J (1) with /j = i if Xy = 
1. Generate F (f ) via formula (10). Define x* = (jq ,...,x* ; ) w *th 

x* =F t (f). 

Step 11: Generate an optimal solution. 
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Combine x* with the solution obtained from (5) and (6) to yield 
an optimal solution of problem (3). 

For more refer [58]. 



1.5 Composite FRE-resolution based on Archimedean triangular 
norms 

The resolution problem of FRE is one of the most important and 
widely studied problems in the field of fuzzy sets and fuzzy 
systems. The first step for the resolution of a FRE is to establish 
the existence of the solution. If the equation is solvable the 
solution set contains a maximum solution and possibly several 
minimum solutions. It has been proved that the finding of these 
solutions suffices for the finding of solution set. 

[92] present sup t FREs. They prove in most practical cases 
the solution set of sup t FREs is non-empty and provide some new 
criteria for checking the existence of the solution. 

Introducing the ‘ solution matrices’ formulation of the 
problem, we find an “if and only if’ condition for the solution 
existence. Then, after a brief description of the most convenient 
algorithm for solving sup-t FREs proposed by [9], a fast algorithm 
is given for determining the solution set of sup-t FREs, where t is 
an Archimedean t-norm. 

Let X, Y, Z be discrete crisp sets with cardinalities n, m and 
k, respectively, and A(X, Y), R(Y, Z), B(X, Z) be three binary 
fuzzy relations constraining its other with the relationship. 

A o' R = B, (1) 

where o' is the well-known sup-t composition (t is a triangular 
norm). Eq. (1) can be written in the matrix form 

A o'R = B, (2) 

where A„ x m , R m x k , B n x k are the matrix representations of A, R, 
B, respectively. Eq. (2) is the typical form of a fuzzy relation 
equation (FRE) for which the following problems arise: 

(i) the resolution of (2) for R, when A and B are 
known, and 
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(ii) the resolution of (2) for A, when R and B are known 
(inverse problem). 

It is remarked that if we have a method for solving the first 
problem, using the same method for the equation R' 1 o' A’ 1 = B' 1 
that employs transposed matrices, the second problem could be 
solved. Thus, without loss generality, we consider only the first 
problem. Moreover, (2) is actually a set of k simpler fuzzy 
relation equations that can be solved independently and so it 
suffices to consider only the equation 

Ao*r = b, (3) 

Where r mx] and b nx i are column vectors of R and B, respectively. 
Clearly, (3) is a system of n equations of the form 

a o* r = b, (4) 

where a is a row vector of A. One can easily see that (3) has 
solution for r if and only if (iff) all the n equations of form (4) 
have at least one common solution for r. 

Let S (A, b) be the solution set of (3), i.e. S(A, b) = (r : 
A o' r = b}. It is well known that if S(A, b) * <|>, then it contains a 
unique maximal solution r and may contain several minimal 
solutions r [43]. The solution set is the union of all the lattices 
[ r , r ] between each minimal and the maximum solution. 




In Figure 1.5.1, a view of the solution set is illustrated. We 
define the mean solution, as the minimal element of the 
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intersection of the lattices [r ,?]. Obviously, the mean solution 
always exists (since r always belongs to the intersection of the 
lattices [ r , r ]) and it is unique. 

The maximum, the mean and the minimum solutions of (2) 
come from the respective solutions of (3) with the aid of the 
following: 

(i) The maximum solution is the m x k matrix 

[ r i ? i...r J, (5) 

where r, (i = 1,2,. . .,k) is the maximum solution of the ith equation 
of form (3). 

(ii) The mean solution is the m x k matrix 

[h ? 2 r k ], (6) 

where r t (i = l,2,...k) is the mean solution of the i th equation of 
form (3). 

(iii) The minimum solutions are the m x k matrices 

[r i r 2 ... r k ], (7) 

where r , e 5, (i = 1,2,..., k) and S i is the minimal solutions set 
of the i th equation of form (3). 

The basis on which equations of form (3) can be solved is the 
simple equation of the form. 

t (a, x) = b (8) 

where a and b are given Eq. (8) is actually a special case of (3), 
for n = m = 1. Our purpose in this section is the study of (8). 

A function t = [0, 1] x [0, 1] — > [0, 1] is a t-norm iff Va, b, de [0, 
1] it satisfies the following four axioms (axiomatic skeleton for t- 
norms): 

Axiom 1: t (a, 1) = a and t (a, 0) =0 
(boundary condition). 
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Axiom 2: b < d implies t (a, b) < t (a, d) 

(monotonicity). 

Axiom 3: t (a, b) = t (b, a) 

(commutativity). 

Axiom 4: t (a, t (b, d)) = t (t (a, b), d) 

(associativity). 

A t-norm t is called Archimedean iff 

Axiom 5: t is a continuous function. 

Axiom 6: t (a, a) < a, Va e (0, 1) 

(subidempotency). 

The class of t-norms has been widely studied by many 
researchers. It is proved in [43] that min (a, b) is the only 
idempotent t-norm. On the other hand, for any t-norm t it is true 
that [43] 



t (a, b) < min (a, b). 

Thus, the only continuous t-norm that is not Archimedean is the 
min (a, b). This is a very interesting result, since it suggests a 
separate study of (3) when t is an Archimedean t-norm and when t 
is the min (a, b). 

Equations of the form (8) do not always have a solution and 
when they do have one, it may not be unique. The following 
proposition can be easily proved. 

Proposition 1.5.1: Let t be a continuous t-norm and a, b, x e [0, 
1]. The equation t(A, x) = b has a solution for x iff a > b. 

Let a® f b and a® f b denote the maximal and the minimal, 
respectively, solution of (8) (if they exist), i.e. 

a <§>' b = sup {x e [0, 1]: t (a, x) = b}, 
a d>' b = inf {x e [0, 1]: t (a, x) = b}. 
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If the solution is unique, it is denoted by a b. Based on the 
above notations, the maximal solution operator (max-SO) cd, and 
the minimal solution operator. (min-SO) cd, are defined as 
follows: 



- , { !, a<b, 

co, (a, b) = < ~ , 

[a ® b, a> b, 


(9) 


\ 0, a < b, 

CD, (a, b) = \ - 

|a ® b, a> b , 


(10) 



when (8) has a unique solution, min-SO and max-SO take the 
form 



( a > b ) = 



1, a < b, 
[a ®‘ b, a> b. 



m,(a, b) 



10, a < b, 
[a ®' b, a>b, 



Max-SO and min-SO are extensions of the a and the <r operators 
defined by Sanchez [84] for the min t-norm. The max-SO 
extension for any t-norm has been proposed by Di Nola et al [17]. 
[9] proposed the min-SO extension. The definition of max-SO 
given in [43] is 

o>, (a, b) = sup {x e [0, 1]: t (a, x) <b}. (11) 

For any a, b e [0, 1], it is cb min ( a,b ) < a and cb min (a, b) < b. 

For any a, b, d e [0, 1] with a < b , it is ro min ( a,b ) < a> min (b,d) . 
For any a, b e [0, 1] and any continuous t-norm t it is 
t (ad), (a, b)) < b. 

For any a, b e [0, 1] and any continuous t-norm t it is co, (a, b) > 
d>, (a, b). 

Flere the resolution of (3) is studied. First, some conditions 
for the existence of the solution are given and then a resolution 
method is described. [9] has originally proposed the method. It is 
formulated in a different way, based on the solution matrices, in 
order to be more comprehensive. Moreover, some new results are 
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given on the form of the solution matrices as well as on the 
solution of existence problem that help to clear up the underlying 
mechanism of the resolution process of FREs. 

Let us first proceed with the solution existence problem. Eq. 
(3) has a solution iff all its equations of form (4) have a common 
solution. The following lemma can be established. 

Let a o' r = b a FRE of form (4). We have S (a, b) ^(J> iff there 
exists j g N m such that aj>b. 

Let A o' r = b be FRE of form (3) with m > n. If for any i g 
N n there exists j g N m such that Ay > b; and A kj - < b k , Vk g N„ - 
{I}, then S (A, b) ^ (|). 

Let t be a continuous t-norm and A o* r = b be a FRE of form 
(3). The matrix T mxn is the mean solution matrix (mean -SM) of 
the FRE, where 

fy = oi, (Aji.bj), Vi e N m ,j g N m 

The matrix E is the maximal solution matrix (max-SM) of the 
FRM, where 

fy = 01, (Ajj bj), Vi G N m ,j G N„, 

The matrix E is the minimal solution matrix (min-SM) of the 
FRM, where 

fy = co, (Ajybj), Vi g N m ,j g N„, 

The matrix f is the minimal solution matrix (min-SM) of the 
FRE, where 

fy = ra min inf f ik , If , Vi g N m , j g N n , 

\ y k^N n J 

Let t be a continuous t-norm and A o' r = b be a FRE of form (3). 
If S (Aj„ bj) * (|i for any j g N„, then 

r. ; eS (Aj.,bj),Vj sN n 
T.j G 5 (. Aj.,bj ), V; G N n 
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where T and F is the min-SM and the mean-SM of the equation, 
respectively. 

Let A o' r = b be a FRE of form (3), where t is a continuous t- 
norm. The following propositions are equivalent: 

(i) S (A, b) = <)>; 

(ii) There exists j e N n , such that A. . =0 and bj 0. 



Algorithm 1 

Step 1: Write down the pseudo-polynomial form of min-SM: 



p= I IE— ( 12 ) 

j£N„i<EN„, 1 
bj* 0 

If r,. ; = 0, it is omitted from (12). All the operations involved in 
(12) (summation, multiplication, division) are symbolic. 

Step 2: Calculate P according to the polynomial multiplication in 
symbolic form. 

Step 3: Simplify P by multiplication and then summation. 

Multiply the terms of the sum using the formula 



a b ( 1TOX ia ' b) l = k, 

T’T = 1 / 

[unchanged, l -t- k. 

Sum among the terms using the formula 

C| Cp c p + dy d 2 dp 

h h ^1 ^2 kp 



£2_ c _p_ 

h ’ h 1 P 

unchanged 



c t <d h \/isN p , 
otherwise. 



(13) 



(14) 
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Step 4: Suppose that after the Step 3, P has s terms. Then (3) has s 
minimal solutions. They are computed using the following 
equation: 

(15) 

where rj° = c ( /K j = 1, 2, ...,m, i= 1, 2,..., s. 

Here some theoretical results are provided that lead to the 
simplification of the method described in the above whenever t is 
an Archimedean t-norm. Note that we have mentioned that the 
only continuous t-norm that is not Archimedean is the 
“minimum”. 

Let us first proceed with some issues on Archimedean 
t-norms. A very important way for generating Archimedean 
t-norms is based on the so-called decreasing generators. A 
decreasing generator is a continuous and strictly decreasing 
function /: [0, 1] — > R such /(l) = 0. 

The pseudo-inverse of/is a function f (-1) : R — » [0, 1] defined 
as 

1, ae (-oo,0), 

/ H) (a) = | f ( ~ r> (a), a e [0,/(0)], (16) 

0, ae (/( 0),+oo) 

where/ -1 is the classical inverse off 

A binary operator t : [0, 1] x [0, 1] — > [0, 1] is an 
Archimedean t-norm iff there exists a decreasing generator /such 
that t(a, b) =/ (-1) (f( a) +/(b)), Va, b e [0, 1], 

Let t be an Archimedean t-norm and a, b, x e [0,1]. The 
equation t(a, x) = b has a solution for x iff a > b. If b e(0, 1) the 
solution x 0 is unique and x 0 e (0,1]. Left as an exercise for the 
reader to prove. 

Algorithm 2 

Step 1 : Write done the pseudo-polynomial form off: 

P= II X j (17) 

yelV„ jsN m 
h ' l ' r r 
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Step 2: Calculate P according to the polynomial multiplication in 
symbolic form. 

Step 3: Simplify P by multiplication and then summation. 

Multiply the sum using the formula 



{l, l =k, 

j/.£, l ^ k 



(18) 



Sum among the terms using the formula 

li- l 2 ...lp + ki.k 2 ...k p 

= | hh-lp if Vi ' G N p ,3jN q : = kj 

[unchanged otherwise. 



Step 4: Suppose that after 3, P has s terms, then (3) has s minimal 
solutions, computed by 



r i j = 1 ’ 

0 otherwise 



j e N m . 



( 20 ) 



We will now show the credibility of the above method. 

The column vector r computed is the maximum solution of 
(3), when t is an Archimedean t-norm. 

The column vector r so found is the mean solution of (3), when t 
is an Archimedean t-norm. 

Algorithm 2 computes the minimal solution set of (3), if t is an 
Archimedean t-norm. 



I. 6 Solving non-linear optimization problem with FRE constraints 

J. Lu, S.C. Fang [61] have used fuzzy relation equation constraints 
to study the non-linear optimization problems. They have 
presented an optimization model with a nonlinear objective 
function subject to a system of fuzzy relation equations. 

The study of the fuzzy relation equations 
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x o A = b 



( 1 ) 



where A = (a ; j ) m x n , 0 < < 1 is a fuzzy matrix lb = (bi,..., 

b n ), 0 < bj < 1 is an n-dimensional vector and ‘0’ is the max-min 
composition [117]. 

The resolution of the equation x o A = b is an interesting and 
on-going research topic. [61] in this paper instead of finding all 
solution of x o A = b, let f(x) be the user’s criterion function, they 
solve the following non linear programming model with fuzzy 
relation constrains 



min f(x) s.t x o A = b (2) 

A minimizer of Eq. (2) will provide a “best” solution to the 
user based on the objective function f(x). some related 
applications of this model with different objective functions can 
be found in [107] for medical diagnosis, and in [60] for 
telecommunication equipment module test. 

Contrary to the traditional optimization problems [62], 
problem (2) subjects to fuzzy relation constraints. From [34], we 
know that when the solution set of the fuzzy relation equations (1) 
is not empty, it is in general a non-convex set that can be 
completely determined by one maximum solution and a finite 
number of minimal solutions. Since the solution set is non- 
convex, conventional optimization methods [60] may not be 
directly employed to solve the problem (2). Recently, [24] studied 
problem (2) with a linear objective function subject to a system of 
FREs and presented a branch and bound procedure to find an 
optimal solution. Flere, we focus on problem (2) with a nonlinear 
objective function and call it a nonlinear optimization problem 
with fuzzy relation constraints (NFRC). A genetic algorithm is 
proposed for solving MFRC. It is designed to be domain specific 
by taking advantage of the structure of the solution set of FREs. 
The individuals from the initial population are chosen from the 
feasible solution set. The genetic operations such as mutation and 
crossover are also kept within the feasible region. It is the beauty 
of this genetic algorithm to keep the search inside of the feasible 
solution set. The well-maintained feasibility of the population 
makes the search more efficient. 

Genetic algorithms (GAs) are built upon the mechanism of 
the natural evolution of genetics. GAs emulate the biological 
evolutionary theory to solve optimization problems. In general. 
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GAs start with a randomly generated population and progress to 
improve solutions by using genetic operators such as crossover 
and mutation operators. In each interaction (generation), based on 
their performance (fitness) and some selection criteria, the 
relatively good solutions are retained and the relatively bad 
solutions are replaced by some newly generated off springs. An 
evaluation criterion (objective) usually guides the selection. 

In the past few years, several methods were proposed to 
handle the constrained optimization problem using genetic 
algorithms. Although there was some variation in details among 
these algorithms, most of them used the penalty or barrier method 
[38, 41]. [64-66] introduced some special genetic operators to 
handle the constrained optimization problem, but those operators 
only work for the problems with a convex domain. A genetic 
algorithm for optimization problem with fuzzy relation constraints 
(GAOFRC) was proposed by them. 

Unlike a general-purpose genetic algorithm, the proposed 
GAOFRC is designed specially for solving nonlinear optimization 
problem with fuzzy relation equations is non-convex in general, 
and the feasible domain is only a small portion of the convex hull 
of it, no existing method is readily available for solving NFRC, 
with the structure inside its feasible domain. 

The proposed GAOFRC uses floating-point representation for 
individuals. Instead of randomly generating a population, the 
initialization process generates a feasible population utilizing the 
structure of fuzzy relation equation. The genetic operators are 
then designed to keep the feasibility of the individuals while they 
evolve. Those solutions with better objective function values will 
have higher opportunities to survive in the procedure. The 
algorithm terminates after it takes a pre-determined number of 
generations. 

In GAOFRC, we use the floating point representation in 
which each gene or variable X; in an individual x = (x b x 2 ,..., x m ) 
is real number from the interval [0 1] since the solution of fuzzy 
relation equations are nonnegative numbers that are less than one. 
More specifically individual x — » (x b ..., x m ) where x ; e [0, 1], 
i= 1,2, . .., m. 

Compared to the GAs which have no feasibility requirement, 
GAOFRC ’s feasibility of individuals limits the search process to a 
much smaller space. 

In general, a GA initializes the population randomly. It works 
well when dealing with unconstrained optimization problems. 
However, for a constrained optimization problem, randomly 
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generated solutions may not be feasible. Since GAOFRC intends 
to keep the solutions (Chromosomes) feasible, we present an 
initialization module to initialize a population by randomly 
generating the individuals inside the feasible domain. 

Since some elements will never play a role in determining the 
solutions to fuzzy relation equations. Therefore, we can modify 
the fuzzy relation matrix by identifying those elements and 
changing their values to us with the hope of easing the procedure 
of fining a new solution. To make it clear, we define some 
“equivalence operators”. 

DEFINITION 1.6.1: If a value-changing iii the element(s) of a 
given fuzzy relation matrix A has no effect on the solutions of 
fuzzy relation equations (1). This value changing is called an 
equivalence operation. 

Lemma 1.6.1: For j t , j 2 e {1, 2,...,n}, if b j> b j ,ay>b ^ and 
a ih > bj for some i, then “resetting a ih to zero” is an 
equivalence operation. 

Based on this idea, the initialization module originates a 
population consisting of a given number of randomly generated 
feasible solutions. The algorithm for initializing is described as 
follows: 









where 



Get the matrix A, b, and size of population P size . 
Compute the potential maximum solution X as follows: 



x = ( A@b ) = 



A i a ij@ b j ) 

7=1 



( 3 ) 



Jl xm 




if ay < bj 
if ay > bj' 



• If x o A = b, continue. Otherwise, stop, the problem is 
infeasible. 

• Simplify matrix A by the equivalence operations. 

• For each element a,j of A, 

Initialize a lower bound parameter lb(i,j) = 0. 

• For i = 1, ...,m, j = 1, ...,n 
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If ay > bj, set lb (i,j) = bj; 

• For i = 1, m, 

Set the maximal lower bound LB max (i) = max"*, lb(i, j). 

• Set k = 1. 

• WHILE (k<P size ) 

Fori= l,...,m, 

Generate a random number pop (k, i) in the interval 

[LFWi), (i)] 

Set k <— k + 1. 

• Output matrix [pop (k,i)] P as the initial population 
of size P s j ze . 

Basically, we took up each individual equation and obtained the 
lower bound for the solutions for each variable X;. Then, we 
compute the maximal lower bound for each x ; . [LB max (i), X (i)] is 
usually an interval. The collection of these intervals we can 
generate a random number in the interval [LB max (i), x (i)] for 
each Xj. 

P(Select the rth individual) = q 1 (l-q) (r_1) , 

where q is the probability of selecting the best individual, r is the 
rank of the individual, q’ = q/(l- (1 - q) Ps,ze ), and P S i ze is the 
population size. 

Because GAOFRC would like to stay feasible, we cannot 
mutate the chromosomes randomly. Although various mutation 
operators handling the constrained optimization problems have 
been proposed in the literature [38, 41, 61, 65], they are all 
designed for the convex problems. There seldom is any mutation 
operator available. 

Table 1 

Elements for each columns such that aij>bj 



Column 1 


Column 2 


Column 3 


Column 4 


Column 5 


a 21, a 31 


a 32 


a 43 


a 14, a 24, a 54, 


a 45 



for the non-convex problem. In what follows, we present a 
mutation operator for GAOFRC whose feasible domain is non- 
convex. 
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Note that a chromosome in GAOFRC is represented by a 1 x 
m vector x = (x [5 x 2 ,..„ x m ). For a given chromosome x 1 = 
[x\,x\,--- ,x ] n ), we define a feasible mutation operator that 
mutates the chromosome by randomly choosing i 0 from 1, 2,...,m 
and decreasing x :} to a random number between [0, xj ], while 

this operation may make the chromosome x 1 infeasible we can 
adjust other x'j. i ^ i 0 , to make x 1 feasible. In fact, the adjustment 
of making the infeasible solution become feasible in nothing but a 

process of finding a new solution. When the changing of xj o pulls 
the x 1 outside of the convex hull of the feasible domain, 
decreasing xj o results in an infeasible solution no matter how 

other x 1 ; s are adjusted. In this case, GAOFRC will neglect this 
decreasing operation and find another x 1 ; to decrease. Since both 
the choosing of x 1 , and the extent of decreasing are randomly 
done, it is guaranteed that a feasible mutation is eventually 
attainable. We present a feasible mutation operation as follows: 

1. Get the simplified matrix A, b and x = (x,, x 2 ,. . ., x m ). 

2. Find the decrease set D, a subset of {1, 2,...,n}, such that 
there are more than one ay at column j of A satisfying 
that ay > bj, for ieD. 

3. Randomly choose an element k from D; generate a 
random number x’ k from the interval [0, x k ], set x <— (xj, 
x 2 ,„, x’ k ,...,x m ). 

m 

4. If (x ; A ay) = bj, Vj, go to step 7; otherwise, go to 

i = 1 

next step. 

5. Generate the increase set N = D - {k}. 

6. For an equation j in which x is not satisfied, randomly 
choose an element x t from the increase set N such that Xj 
< bj and ay > bj. Set Xj = bj, go to step 4. 

7. Go to crossover operation. 
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